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LAUNCHING OF POYNTING JETS FROM ACCRETION DISKS 

R. V. E. Lovelace, 1 M. M. Romanova, 1 
RESUMEN 

El resumen sera traducido al espanol por los editores. The jets observed to emanate from many 
compact accreting objects may arise from the twisting of the magnetic field threading a differentially rotating 
accretion disk which acts to magnetically extract angular momentum and energy from the disk. Two main 
regimes have been discussed, hydromagnetic outflows, which have a significant mass flux and have energy and 
angular momentum carried by both matter and electromagnetic field and, Poynting outflows, where the mass 
flux is negligible and energy and angular momentum arc carried predominantly by the electromagnetic field. We 
describe recent theoretical work on the formation of relativistic Poynting jets from magnetized accretion disks 
and new relativistic, fully-electromagnetic, particle-in-ccll simulations of the formation of jets from accretion 
disks. 

ABSTRACT 

The jets observed to emanate from many compact accreting objects may arise from the twisting of the magnetic 
field threading a differentially rotating accretion disk which acts to magnetically extract angular momentum 
and energy from the disk. Two main regimes have been discussed, hydromagnetic outflows, which have a 
significant mass flux and have energy and angular momentum carried by both matter and electromagnetic 
field and, Poynting outflows, where the mass flux is negligible and energy and angular momentum are carried 
predominantly by the electromagnetic field. Wc describe recent theoretical work on the formation of relativis- 
tic Poynting jets from magnetized accretion disks and new relativistic, fully-clcctromagnctic, particle-in-cell 
simulations of the formation of jets from accretion disks. 
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1. GENERAL 

Powerful, highly-collimated, oppositely directed 
jets are observed in active galaxies and quasars, 
and in old compact stars in binaries - the 'micro- 
quasars'. Different models have been put forward 
to explain astrophysical jets (Bisnovatyi-Kogan & 
Lovelace 2001). Recent observational and theoret- 
ical work favors models where twisting of an ordered 
magnetic field threading an accretion disk acts to 
magnetically accelerate the jets. Two main regimes 
have been considered in theoretical models, the hy- 
dromagnetic regime where the energy and angular 
momentum is carried by both the electromagnetic 
field and the kinetic flux of matter, and the Poynt- 
ing flux regime where the energy and angular mo- 
mentum outflow from the disk is carried predomi- 
nantly by the electromagnetic field. In §3 we out- 
line the theory of Poynting jets. In §4 we present 
new results from axisymmctric, fully electromag- 
netic, relativistic-particle-in-cell (PIC) simulations 
of the formation and propagation of relativistic jets 
from a disk. 



1 Cornell University, Ithaca, NY USA. 



2. THEORY OF POYNTING JETS 

The powerful jets observed from active galaxies 
and quasars are probably not hydromagnetic out- 
flows but rather Poynting flux dominated jets. The 
motion of these jets measured by very long baseline 
interferometry correspond to bulk Lorentz factors of 
F = O(10) which is much larger than the Lorentz 
factor of the Keplcrian disk velocity predicted for 
hydromagnetic outflows. Furthermore, the low Fara- 
day rotation measures observed for these jets at dis- 
tances < kpc from the central object implies a very 
low plasma densities. Similar arguments indicate 
that the jets of microquasars are not hydromagnetic 
outflows but rather Poynting jets. Poynting Jets 
have been proposed to be the driving mechanism 
for gamma ray burst sources (Katz 1997). Theo- 
retical studies have developed models for Poynting 
jets from accretion disks (Lovelace, Wang, & Sulka- 
nen 1987; Lyndcn-Bcll 2003; Romanova & Lovelace 
1997; Levinson 1998; and Lovelace et al. 2002; and 
Lovelace & Romanova 2003 ). Stationary Poynting 
flux dominated outflows were found by Romanova et 
al. (1998) and Ustyugova et al. (2000) in axisym- 
metric MHD simulations of the opening of magnetic 
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Fig. 1. Sketch of the magnetic field configuration of a 
Poynting jet from Lovelace and Romanova (2003). The 
bottom part of the figure shows the initial dipole-like 
magnetic field threading the disk which rotates at the 
angular rate fi(r). The top part of the figure shows the 
jet at some time later when the head of the jet is at 
a distance Z(t). At the head of the jet there is force 
balance between electromagnetic stress of the jet and the 
ram pressure of the ambient medium of density p ex t- 

loops threading a Keplerian disk. 

We first summarize the theory of non-relativistic 
Poynting jets which is based on the Grad-Shafranov 
equation. We show results of non-relativistic, ax- 
isymmetric MHD simulations which support this 
theory. Later, we discuss the corresponding results 
obtained by solving the relativistic Grad-Shafranov 
equation. 

Consider the coronal magnetic field - such as that 
shown in the lower part of Figure 1 - of a differ- 
entially rotating Keplerian accretion disk. That is, 
the disk is perfectly conducting, high-density, and 
has a small accretion speed (<C dr-). Further, con- 
sider "coronal" or "force-free" magnetic fields in the 



non-relativistic limit. We use cylindrical (r, <f>, z) 
coordinates and consider axisymmetric field config- 
urations. Thus the magnetic field has the form 
B = Bp + B^cj) , with B p = B r r + B z z. Be- 
cause V B = 0, B = VxA with A the vector 
potential. Consequently, B r = — (l/r)d^/dz and 
B z = {l/r)d^/dr. where *(r, z) = rA<j,(r,z). The 
fy(r, z) = const lines label the poloidal field lines; 
that is, (Bp • V)* = = (B • V)*. Note that 
(r, z) is the magnetic flux through a horizontal, 
coaxial circular disk of radius r. The magnetic field 
threading the disk at z = is assumed to evolve 
slowly so that it can be considered approximately 
time- independent, ^(r, z = 0) = ^q(t). However, 
the magnetic field above the disk will in general be 
time-dependent, <!> = ty(r,z,i), due to the differen- 
tial rotation of the disk. 

The non-relativistic equation of plasma motion 
in the corona of an accretion disk is pdw/dt = 
—\7p+pg + J x B/c, where v is the flow velocity, p is 
the pressure, and g is the gravitational acceleration. 
The equation for the B field is V x B = 4ttJ/c, be- 
cause the displacement current is negligible in the 
non-relativistic limit. In the coronal or force- free 
plasma limit, the magnetic energy density B 2 /87r 
is much larger than the kinetic or thermal energy 
densities; that is, we have sub-Alfvenic flow speeds 
v2 ^ V A — B 2 /(4 7r /o), where va is the Alfven veloc- 
ity. The force equation then simplifies to w J x B 
so that J = AB (Gold & Hoylc 1960). Because 

V • J = 0, (B ■ V)A = 0, and consequently A = 
A(\t), as well-known. Thus Ampere's law becomes 

V x B = 47rA(^)B/c. The r and z components of 
this equation imply rB^ = H(ty), and dH{^)/d^> = 
47rA( , F)/c, where is another function of 
Thus, H(^>) = const are lines of constant poloidal 
current density; J p = (c/4ir)(dH/d^f)'Bp so that 
(Jp • W)H = 0. The toroidal component of Ampere s 
law gives the non-relativistic Grad-Shafranov equa- 
tion for 

A** = . (i) 

Here, A* = d 2 /dr 2 - (l/r)(d/dr) + d 2 /dz 2 is the ad- 
joint Laplacian operator. Note that A*^ = r(V 2 — 
l/r 2 )A and that H(dH/dV) = AnrJ^/c. From 
Ampere's law, § d\ ■ B = (Att/c) J dS ■ J, so that 
rB,p(r,z) = H(^) is (2/c) times the current flowing 
through a circular area of radius r (with normal z) la- 
beled by $(r, z) = const. Equivalently, — H [^(r, 0)] 
is (2/c) x the current flowing into the area of the disk 
with radii < r. For all cases studied here, —H(^f) 
has a maximum so that the total current flowing into 
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the disk for r < r m is I = (2/c)(-H) max , where r m 
is such that — H[^(r m , 0)] = (—H) max so that r m is 
less than the radius of the O-point, ro. The same 
total current hot flows out of the region of the disk 
r = r m to r . 

The function H(ty) must be determined before 
the Grad-Shafranov equation can be solved. H(^) 
is determined by the differential rotation of the disk: 
The azimuthal twist of a given field line going from 
an inner footpoint at r\ to an outer footpoint at ri 
is fixed by the differential rotation of the disk. The 
field line slippage speed through the disk due to the 
disk's finite magnetic diffusivity is estimated to be 
negligible compared with the Keplerian velocity vk ■ 
For a given field line we have rdc^/B^ = ds p /B p , 
where ds p = \J dr 2 + dz 2 is the poloidal arc length 
along the field line, and B p = B 2 . + B 2 . The total 
twist of a field line loop is 




with the sign included to give A<fi > 0. For a Keple- 
rian disk around an object of mass M, the angular 
rotation rate is SIk = \J GM/r 3 so that the field 
line twist after a time t is A</>(<3>) = Q t [(r /ri) 3 / 2 — 
(r /r 2 ) 3/2 ] = (Slot)F(V/V ), where r is the radius 
of the O-point, ft = y/GM/rl, and F is a dimen- 
sionless function. 

Equations (1) and (2) have been solved numeri- 
cally by Li et al. (2001) and Lovelace et al. (2002) 
for an initial poloidal magnetic field as shown in the 
lower part of Figure 1. As the "twist," as measured 
by Qot, increases, a high twist field configuration ap- 
pears with a different topology. A "plasmoid" con- 
sisting of toroidal flux detaches from the disk and 
propagates outward. The plasmoid is bounded by a 
poloidal field line which has an X-point above the 
O-point on the disk. The occurrence of the X-point 
requires that there be at least a small amount of 
dissipation in the evolution from the poloidal dipolc 
field and the Poynting jet configuration. The high- 
twist configuration consists of a region near the axis 
which is magnetically collimated by the toroidal B^ 
field and a region far from the axis which is anti- 
collimatcd in the sense that it is pushed away from 
the axis. The field lines returning to the disk at 
r > ro are anti-collimated by the pressure of the 
toroidal magnetic field. The poloidal field fills only 
a small fraction of the coronal space. 

Figure 2 shows results of non-relativistic axisym- 
metric simulations of the formation of a Poynting 
jet by Ustyugova et al. (2000). The flow near the 
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Fig. 2. Non-relativistic time evolution of dipole-like 
field threading the disk from the initial configuration 
t = (bottom panels) to the final quasi-stationary state 
at t = 1.2t out , where t out is the rotation period of the 
disk at the outer radius R ou t of the simulation region 
from Ustyugova et al. (2000). The left hand panels 
show the poloidal field lines which are the same as the 
ty(r,z) =const lines; ^ is normalized by ^ m ax, and the 
spacing between lines is 0.1. The middle panels show the 
poloidal velocity vectors v p . The right-hand panels show 
the constant lines of —rB^(r, z) > in units of $ mal /ro, 
and the spacing between lines is 0.1. For this calcula- 
tion a 100 x 100 inhomogeneous grid was used with Arj 
and Azfc growing with distance r and z geometrically 
as Arj = Ariq' and Azk = Az\q k , with q = 1.03 and 
An = Azi = 0.05r . 

z— axis is the Poynting jet and its physical properties 
agree with Grad-Shafranov solutions. 

In the case of relativistic Poynting jets we hy- 
pothesize that the magnetic field configuration is 
similar to that in the non-relativistic limit (Ustyu- 
gova et al. 2000; Lovelace et al. 2002). Thus, 
most of the twist A<fi of a field line of the relativis- 
tic Poynting jet occurs along the jet from z = to 
Z(t) as sketched in Figure 3, where Z(t) is the ax- 
ial location of the "head" of the jet. Along most of 
the distance z = to Z : the radius of the jet is a 
constant and 4* = 4"(r) for Z >> ro. Note that 
the function ^f(r) is different from ^(r, 0) which is 
the flux function profile on the disk surface. Hence 
r 2 d(j)/dz = rB^(r,z)/B x (r,z). We take for simplic- 
ity V z = dZ/dt = const. We determine V z subse- 
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Fig. 3. Sketch of the initial conditions for the relativis- 
tic particle-in-cell simulations of jet formation from an 
accretion disk. 

quently. In this case H{9) = [r 2 tt(^)/V z ]B z can 
be written as a function of VP and d^Sf/dr. With H 
known, the relativistic Grad-Shafranov equation, 
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r 4 ft 2 



dH{^) 
d^ ' 
(3) 

can be solved (Lovelace & Romanova 2003). 

The quantity not determined by equation (3) 
is the velocity V z , or Lorentz factor T = 1/(1 — 
V z / c 2 ) 1 / 2 . This is determined by taking into ac- 
count the balance of axial forces at the head of the 
jet: the electromagnetic pressure within the jet is 
balanced against the dynamic pressure of the exter- 
nal medium which is assumed uniform with density 
Pext . This gives (r 2 - l) 3 = B 2 /(8TrK 2 Pext c 2 ), or for 
F> 1, 



r 



1/3 



/_Bo_y /3 

\10 3 GJ 



l/cm a 

Hext 



1/6 



(4) 



where 7Z = ro/r g ^> 1, with ro the O-point of the 
magnetic field, r g = GM/c , and Bq the magnetic 
field strength at the center of the disk. This value of 
r is of the order of the Lorentz factors of the expan- 
sion of parsec-scale extragalactic radio jets observed 
with very-long-baseline-interferometry (see, e.g., 
Zensus et al. 1998). This interpretation assumes 



that the radiating electrons (and/or positrons) are 
accelerated to high Lorentz factors (7 ~ 10 3 ) at the 
jet front and move with a bulk Lorentz factor F rela- 
tive to the observer. The luminosity of the +z Poynt- 
ing jet is Ej = c /J" 2 rdrE r B^I1 = cB 2 U 3 / 2 r 2 /3 ~ 
2.1 x 10 46 (£ /10 3 G) 2 (ft/10) 3 / 2 (M/10 9 M Q ) 2 erg/s, 
where M is the mass of the black hole. 

For long time-scales, the Poynting jet is of course 
time-dependent due to the angular momentum it ex- 
tracts from the inner disk (r < ro) which in turn 
causes ro to decrease with time (Lovelace et al. 
2002). This loss of angular momentum leads to a 
"global magnetic instability" and collapse of the in- 
ner disk (Lovelace et al. 1994, 1997, 2002) and a 
corresponding outburst of energy in the jets from 
the two sides of the disk. Such outbursts may ex- 
plain the flares of active galactic nuclei blazar sources 
(Romanova & Lovelace 1997; Levinson 1998) and 
the one-time outbursts of gamma ray burst sources 
(Katz 1997). 

3. RELATIVISTIC PARTICLE-IN-CELL 
SIMULATIONS OF JETS 

Wc performed relativistic, fully electromagnetic, 
particle-in-cell simulations of the formation of rela- 
tivistic jets from an accretion disk initially threaded 
by a dipolc-likc magnetic field. This was done using 
the code XOOPIC developed by Verboncoeur, Lang- 
don, and Gladd (1995). Earlier, Gislcr, Lovelace, 
and Norman (1989) studied jet formation for a 
monopole type field using the relativistic PIC code 
ISIS. The geometry of the initial configuration is 
shown in Figure 3. The computational region is a 
cylindrical "can," r = — R m and z = — Z m , with 
outflow boundary conditions on the outer bound- 
aries, and the potential and particle emission speci- 



fied on the disk surface r = — R„ 



0. Equal 



fluxes of electrons and positrons are emitted so that 
the net emission is effectively space-charge-limitcd. 
About 10 5 particles were used in the simulations re- 
ported here. The behavior of the lower half-space 
(z < 0) is expected to be a mirror image of the up- 
per half-space. 

Figure 4 shows the formation of a relativistic jet. 
The gray scale indicates the logarithm of the den- 
sity of electrons or positrons with 20 levels between 
the lightest (10 12 ) and darkest (4 x 10 15 /ru 3 ). The 
lines are poloidal magnetic field lines B p . The to- 
tal, three-dimensional magnetic field is shown in Fig- 
ure 5. The computational region has (R m , Z m ) — 
(50, 100) m, the initial B— field is dipole-like with 
5,(0,0) = B = 28.3 G and an O-point at {r, z) = 
(10,0) m, and the electric potential at the center of 
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Fig. 4. Relativistic particle-in-cell simulations of the 
formation of a jet from a rotating disk, (a) -(c) give 
snapshots at times (1, 2, 3) x 10~ 7 s, and (d) is at t = 1CP 6 
s. 

the disk is <t>o = — 10 7 V relative to the outer region 
of the disk. Initially, the computational region was 
filled with a distribution of equal densities of elec- 
trons and positrons with n±(0, 0) = 3 x 10 13 /m 3 . 
Electrons and positrons are emitted with equal cur- 
rents I± = 3 X 10 5 A from both the inner and the 
outer portions of the disk as indicated in Figure 3 
with an axial speed much less than c. For a Kcplcrian 
disk with r S> r g , the scalings are 4> ~ B^r^rg) 1 ' 2 , 
I ~ cBqTq and the jet power is ~ cB^r^ 2 ^ 2 . The 
calculations were done on a 64 x 128 grid stretched in 
both the r and z directions so as to give much higher 
spatial resolution at small r and small z. These sim- 
ulations show the formation of a quasi-stationary, 
collimatcd current-carrying jet. The Poynting flux 
power of the jet is Ej ~ 7 x 10 11 W and the particle 
kinetic energy power is r* 4.7 x 10 10 W. The charge 
density of the electron flow is partially neutralized 
by the positron flow. Simulations arc planned with 
the positrons replaced by ions. 

We thank the meeting organizers for the stimu- 
lating and very well organized meeting. This work 
was supported in part by DOE cooperative agree- 
ment DE-FC03 02NA00057. 
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